Deflection of field-free aligned molecules 
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We consider deflection of polarizable molecules by inhomogeneous optical fields, and analyze the 
role of molecular orientation and rotation in the scattering process. It is shown that molecular rota- 
tion induces spectacular rainbow-like features in the distribution of the scattering angle. Moreover, 
by pre-shaping molecular angular distribution with the help of short and strong femtosecond laser 
pulses, one may efficiently control the scattering process, manipulate the average deflection angle 
and its distribution, and reduce substantially the angular dispersion of the deflected molecules. This 
opens new ways for many applications involving molecular focusing, guiding and trapping by optical 
and static fields. 

PACS numbers: 33.80.-b, 37.10.Vz, 42.65.Re, 37.20.+j 



Optical deflection of molecules by means of nonreso- 
nant laser fields is a hot subject of many recent exper- 
imental studies [H, 0, S 0j- By controlling molecular 
translational degrees of freedom ||, novel elements of 
molecular optics can be realized, including molecular lens 

0, 0] and molecular prism Q . The mechanism of molec- 
ular deflection by a nonuniform laser field is rather clear: 
the field induces molecular polarization, interacts with 
it, and deflects the molecules along the intensity gradi- 
ent. As most molecules have anisotropic polarizability, 
the deflecting force depends on the molecular orientation 
with respect to the deflecting field. Previous studies on 
optical molecular deflection have mostly considered ran- 
domly oriented molecules, for which the deflection angle 
is somehow dispersed around the mean value determined 
by the orientation-averaged polarizability. The latter be- 
comes intensity-dependent for strong enough fields due 
to the field-induced modification of the molecular angu- 
lar motion [f| 0] • This adds a new ingredient for control- 
ling molecular trajectories [1, H|, which is important, but 
somehow limited because of using the same fields for the 
deflection process and orientation control. 

In this Letter, we show that the deflection process can 
be significantly affected and controlled by pre-shaping 
molecular angular distribution before the molecules en- 
ter the interaction zone. This can be done with the 
help of numerous recent techniques for laser molecular 
alignment, which use single or multiple short laser pulses 
(transform-limited, or shaped) to align molecular axes 
along certain directions. Short laser pulses excite rota- 
tional wavepackets, which results in a considerable tran- 
sient molecular alignment after the laser pulse is over, 

1. e. at field-free conditions (for recent reviews on field- 
free alignment, see, e.g. (HQ). Field-free alignment was 
observed both for small diatomic molecules as well as for 
more complex molecules, for which full three-dimensional 
control was realized [TfJ . [Til [T3 | . We demonstrate that 
the average scattering angle of deflected molecules and 
its distribution may be dramatically modified by a proper 
field- free pre-alignment. By separating the processes of 
the angular shaping and actual deflection, one gets a flex- 
ible tool for tailoring molecular motion in inhomogeneous 
optical and static fields. 



Although our arguments are rather general, we fol- 
low for certainty a deflection scheme that reminds the 
experiment by Stapelfeldt et al [l[ who used a strong 
IR laser to deflect a CS2 molecular beam, and then ad- 
dressed a portion of the deflected molecules (at a pre- 
selected place and time) by an additional short and 
narrow ionizing pulse. Consider deflection (in z direc- 
tion) of a linear molecule moving in x direction with 
velocity v x and interacting with a focused nonresonant 
laser beam that propagates along the y axis. The spa- 
tial profile of the laser electric field in the xz-plane is 
E = E a cxp[-(x 2 + z 2 )/uj 2 } exp[-21n2t 2 /r 2 ]. The inter- 
action potential of a linear molecule in the laser field is 
given by: 



U(t) = —E 
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where E is defined above, and a\\ and a± are the com- 
ponents of the molecular polarizability along the molec- 
ular axis, and perpendicular to it, respectively. Here 8 is 
the angle between the electric field polarization direction 
(along the laboratory z axis) and the molecular axis. A 
molecule initially moving along the x direction will ac- 
quire a velocity component v z along z-direction. We con- 
sider the perturbation regime (weak field approximation) 
corresponding to a small deflection angle, 7 « v z /v x . We 
substitute x = v x t, and consider z as a fixed impact pa- 
rameter. The deflection velocity is given by: 
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Here M is the mass of the molecules, and F z is the 
deflecting force. The time-dependence of the force F z (t) 
(and potential U(t)) in Eq.Q comes from three sources: 
pulse envelope, projectile motion of the molecule through 
the laser focal area, and time variation of the angle 
due to molecular rotation. For simplicity, we as- 
sume that the deflecting field does not affect significantly 
the rotational motion. Such approximation is justified, 
say for CS2 molecules with the rotational temperature 
T = 5K, which are subject to the deflecting field of 
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3 • 10 9 W/cm 2 . The corresponding alignment potential 
U « — j (aii — ~ 0-04 me]/ is an order of mag- 

nitude smaller than the thermal energy fc^T, where fee 
is Boltzmann's constant. This assumption is even more 
valid if the molecules were additionally subject to the 
aligning pulses prior to deflection. 

Since the rotational time scale is the shortest one in 
the problem, we average the force over the fast rotation, 
and arrive at the following expression for the deflection 
angle, j = v z /v x : 



7 = 7o [a\\A + a±{l-A)] /< 



(3) 



Here a — l/3ori| + 2/3aj_ is the orientation-averaged 
molecular polarizability, and A = cos 2 9 denotes the 
time-averaged value of cos 2 9. This quantity depends on 
the relative orientation of the vector of angular momen- 
tum and the polarization of the deflecting field. It is 
different for different molecules of the incident ensem- 
ble, which leads to the randomization of the deflection 
process. The constant 70 presents the average deflection 
angle for an isotropic molecular ensemble: 



To 



aE 2 



-4z 



1 



2w 2 In 2 



-1/2 



exp 



2Z 2 
,2 



(4) 



We provide below some heuristic classical arguments on 
the anticipated statistical properties of A and 7 (both 
for thermal and pre-aligned molecules) and then support 
them by a more refined quantum treatment. 

Consider a linear molecule that rotates freely in a 
plane that is perpendicular to the vector J of the 
angular momentum (see Fig.jJTJ). 




FIG. 1: A molecule rotates with a given angular momentum 
J that is randomly oriented in space. 9j is the angle between 
the angular momentum and the laboratory z axis. 

The projection of the molecular axis on the vertical 
^-direction is given by cos 9(t) = cos(wi) sin#j, where 9j 
is the angle between J and z-axis, and 10 is the angular 
frequency of molecular rotation. Averaging over time, 
one arrives at: 
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A = cos 2 9 = — sin 2 1 
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(5) 



In a thermal ensemble, vector J is randomly oriented in 
space, with isotropic angular distribution 1/2 shx(8j)d8j. 



The mean value of the deflection angle is (7) = 70. 
Eq.© allows us to obtain the distribution function, f(A) 
for A (and the related deflection angle) from the known 
isotropic distribution for 9j. Since the inverse function 
9 j (A) is multivalued, one obtains 
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where we summed over the two branches of 9j(A). This 
formula predicts an unimodal rainbow singularity in the 
distribution of the scattering angles at the maximal value 
7 = 7o(ai| + a±)/2a (for A = 1/2), and a flat step near 
the minimal one 7 = 7oa_i_/a (for A = 0). Assume 
now that the molecules are pre-aligned before entering 
the deflection zone by a strong and short laser pulse that 
is polarized perpendicular to the polarization direction 
of the deflecting field (e.g., in re-direction) . Such a pulse 
forces the molecules to rotate preferentially in the planes 
containing the rc-axis. As a result, the vector J of the 
angular momentum is confined to the yz-plane, and an- 
gle 9 j becomes uniformly distributed in the interval [0, n] 
with probability density d9j /n. The corresponding prob- 
ability distribution for A takes the form 
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In contrast to Eq.©, formula Eq.© suggests a bimodal 
rainbow in the distribution of deflection angles, with sin- 
gularities both at the minimal and the maximal angles. 
Finally, we proceed to the most interesting case when the 
molecules are pre-aligned by a short strong laser pulse 
that is polarized parallel to the direction of the deflect- 
ing field. After excitation by such a pulse, the vector of 
the angular momentum of the molecules is preferentially 
confined to the xy-plane, and angle 9j takes a well de- 
fined value of 9j « 7r/2. As a result, the dispersion of 
the scattering angles is reduced dramatically. The distri- 
bution of the deflection angle 7 transforms to a narrow 
peak (asymptotically - a (5-function) near the maximal 
value, 7 = 70 (an + a_i_)/2a. 

For a more quantitative treatment, involving analysis 
of the relative role of the quantum and thermal effects on 
one hand, and the strength of the pre-aligning pulses on 
the other hand, we consider quantum-mechanically the 
deflection of a linear molecule described by the Hamil- 
tonian Ti = J 2 /(2I). Here J is operator of angular mo- 
mentum, and / is the moment of inertia, which is related 
to the molecular rotational constant, B — h/{AirIc) (c is 
speed of light). Assuming again that the deflecting field 
is too weak to modify molecular alignment, we consider 
scattering in different | J, m) states independently. The 
deflection angle is given by Eq. ([3]) , in which A is replaced 
by 



Aj t m — (J, tu\ cos 2 9\J, m) = — + 



2 J(J+ 1) - 3m 2 



3 (2J+3)(2J- 1)' 



(8) 
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In the quantum case, the continuous distribution of the 
angles 7 is replaced by a set of discrete lines, each of 
them weighted by the population of the state | J, to) . 
Fig. [2] shows the distribution of Aj : m in the thermal 
case for various values of the dimensionless parameter 
Jt = \JksT I (hBc) that represents the typical "ther- 
mal" value of J (for Jt > 1). 
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FIG. 2: Quantum distribution of Aj, m in the thermal case. 
Panels (a) and (b) correspond to Jt = 5 and Jt = 15, re- 
spectively. Histogram in panel (c) presents a coarse-grained 
continuous normalized distribution of Aj, m produced from 
(b) by averaging over a set of finite bins. 

The distribution of discrete values of Aj im demon- 
strates a non-trivial pattern, however it shows the ex- 
pected unimodal rainbow feature (see Eq.([6])) for large 
enough Jt after the coarse-grained averaging . 

If the molecules are subject to a strong femtosecond 
prc-aligning pulse, the corresponding interaction poten- 
tial is given by Eq.((T]), in which E(t) is replaced by the 
envelope e(t) of the femtosecond pulse. If the pulse is 
short compared to the typical periods of molecular ro- 
tation, it may be considered as a delta-pulse. In the 
impulsive approximation, one obtains the following rela- 
tion between the angular wavefunction before and after 
the pulse applied at t = 0: 



Hf(t = 0+) = exp(iPcos 2 6»)*(i = 0~ 



(9) 



where the kick strength, P is given by P = (1/AK) ■ (an — 
a ±) J^oo £ 2 (t)dt. Here we assumed the vertical polariza- 
tion (along z-axis) of the pulse. Physically, the dimen- 
sionless kick strength, P equals to the typical amount of 
angular momentum (in the units of K) supplied by the 
pulse to the molecule. For the vertical polarization of 
the laser field, to is a conserved quantum number. This 
allows us to consider the excitation of the states with 



different initial to values separately. In order to find 
fy(t = + ) for any initial state, we introduce an artifi- 
cial parameter £ that will be assigned the value £ = 1 at 
the end of the calculations, and define 

^ = cxp [{iP cos 2 0)£] *(* = 0") =^cj(£)|J,to). 

(10) 

By differentiating both sides of Eq. (fTU|) with respect to 
£, we obtain the following set of differential equations for 
the coefficients cj: 



cj> = cj(J', to| cos 2 6>| J, to), 



(11) 



where c = dc/d^. The diagonal matrix elements in 
Eq. tfTTj) are given by Eq.©, the off-diagonal ones can 
be found using recurrence relations for the spherical har- 
monics [IU- Since * ?=0 = = 0~) and = $>(t = 
+ ) (sec Eq. (fTU|) ). we solve numerically this set of equa- 
tions from £ = to £ = 1, and find #(t = 0+). In 
order to consider the effect of the field-free alignment at 
thermal conditions, we repeated this procedure for every 
initial | Jo, Too) state. To find the modified population of 
the I J, to) states, the corresponding contributions from 
different initial states were summed together weighted 
with the Boltzmann's statistical factors. For symmet- 
ric molecules, statistical spin factor should be taken into 
account. For example, for CS% molecules in the ground 
electronic and vibrational state, only even J values are al- 
lowed due to the permutation symmetry for the exchange 
of two Bosonic Sulfur atoms (that have nuclear spin 0). 




FIG. 3: Distribution of Aj,m for molecules pre-aligned with 
the help of a short laser pulse polarized in the x direction. The 
left column (a-b) presents directly the Aj, m values, while the 
right column (c-d) shows the corresponding coarse-grained 
histograms (as in Fig. [5Ji). Panels (a) and (c) are calculated 
for Jt — 5 and P = 5; (b) and (d) are for Jt = 5 and P — 25. 

Using this technique, we considered deflection of initially 
thermal molecules that were pre-aligned with the help of 
short pulses polarized in x and z directions (Figs. [3] and 
31 respectively) . In the case of the alignment perpendicu- 
lar to the deflecting field, the coarse-grained distribution 
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FIG. 4: Distribution of Aj, m for molecules pre-aligned in the 
z direction. The left column (a-b) presents directly the Aj, m 
values, while the right column (c-d) shows the corresponding 
coarse-grained histograms. Panels (a) and (c) are calculated 
for Jt = 5 and P — 5; (b) and (d) are for Jt = 5 and P = 25. 

of Aj.m (and that of the deflection angles) exhibits the 
bimodal rainbow shape, Eq.0 for strong enough kicks 
(P» 1 and P 3> Jt)- Finally, and most importantly, 
pre-alignment in the direction parallel to the deflecting 
field allows for almost complete removal of the rotational 
broadening. A considerable narrowing of the distribution 
can be seen when comparing Fig. [2^ and Figs. [4)d and 

m. 

Our results indicate that pre-alignment provides an 
effective tool for controlling the deflection of rotating 
molecules, and it may be used for increasing the bright- 
ness of the scattered molecular beam. This might be im- 



portant for nano-fabrication schemes based on the molec- 
ular optics approach Moreover, molecular deflection 
by non-resonant optical dipole force is considered as a 
promising route to separation of molecular mixtures (for 
a recent review, see flij). Narrowing the distribution 
of the scattering angles may substantially increase the 
efficiency of separation of multi-component beams, es- 
pecially when the pre-alignment is applied selectively to 
certain molecular species, such as isotopes [IBf . or nu- 
clear spin isomers [3, E3] • More complicated techniques 
for pre-shaping the molecular angular distribution may 
be considered, such as confining molecular rotation to a 
certain plane by using the " optical molecular centrifuge" 
app roach [l8[ , double-pulse ignited " molecular propeller" 
[19| , or two-direction alignment alternation excited by el- 
liptic laser pulses [20(. In this case, a narrow angular 
peak is expected in molecular scattering, whose position 
is controllable by inclination of the plane of rotation with 
respect to the deflecting field. Laser pre-alignment may 
be used to manipulate molecular deflection by inhomoge- 
neous static fields as well (for recent exciting experiments 
on post-alignment of molecules scattered by static elec- 
tric fields see [HI). In particular, one may affect molec- 
ular motion in relatively weak fields that are insufficient 
to modify rotational states by themselves. Moreover, 
the same mechanisms may prove efficient for controlling 
inelastic molecular scattering off metalic/dielectric sur- 
faces. These and other aspects of the present problem 
are subjects of an ongoing investigation. 
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